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The color connection of the parton system is an important bridge to connect the perturbative 
process and the hadronization one. The special color connection of four-heavy-quark system in e^'e” 
annihilation, which is the necessary one for the doubly heavy baryon and tetraquark productions is 
revisited. The hadronization effects, investigated with the help of the Quark Combination Model are 
compared with the corresponding results employing the Lund String Model in our previous work. 

The global properties related to a certain color connection are not sensitive to various hadronization 
models. 
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In high energy collisions, hadronization is one of the 
most important processes for understanding the non- 
perturbative quantum chromodynamics and the confine¬ 
ment mechanism. e+e“ annihilation has the advantage 
that the hadronization results can be compared with 
the experimental data directly to extract informations of 
hadronization because of no hadrons in the initial states. 
To embed the parton system into a hadronization model, 
it is necessary to specify the concrete color connection. 
However, the color connection for multi-parton system is 
not unique. Some clues can be obtained by analyzing the 
decomposition of the color space of the final partons [ll- 
Q . One of the interesting and important examples is the 
four-heavy-quark system {QQ'QQ'), in which many phe¬ 
nomena related to QCD properties can be studied, e.g., 
the (re)combination of quarks in the production of spe¬ 
cial hadrons, the influence of soft interaction on the re¬ 
construction of intermediate particles (such as W^), etc., 
most of which are more or less related to the color connec¬ 
tions among these four quarks. In Ref. [3, the interesting 
decompositions of the color space of the QQ'QQ' system 
in e+e“ annihilation have been discussed, which are 

(3q 0 3q/) 0 Tq, = 3*qq, <S> 3*q <S> 3*q, ® ■ ■ ■ , or 
3q ® 3q/ (g) (33 (g) 33 ,) = 3q 0 3q/ (g) 3qq, © • • • , (1) 

where 3 and 3* denote the triplet and anti-triplet rep¬ 
resentations of the SUc{3) Group respectively, and the 
subscripts correspond to the relevant (anti)quark. We 
investigated the four-heavy-quark system in which only 
one QQ' (or QQ') pair has a small invariant mass, while 
the invariant masses of the other two are left unre¬ 
stricted, i.e., a three-jet event configuration. When two 
quarks (antiquarks) with small invariant masses in the 
color state 3* (3) attract each other and can be consid¬ 
ered as a diquark (or an antidiquark), it can hadronize 
into a(n) (anti)baryon (tetraquark). The corresponding 
hadronization procedure is a ‘branching’ process via the 
creation of quarks from vacuum by the strong interac¬ 
tions within the system. Let’s take cccc as an example. 


considering the concrete case that (cc) combines with a 
quark q (antidiquark qq') to form Sec [1, 0| {Tec [I])- 
To balance the quantum numbers of color and flavor, an 
antiquark/diquark must be simultaneously created from 
vacuum. To branch them further, more quark pairs and 
diquark pairs must be created from vacuum via the inter¬ 
actions among the quark system. Such a cascade process 
will proceed until the end of time when most of the ‘inner 
energy’ of the entire system is transformed into the kine- 
matical energies and masses of the produced hadrons. 
Each of two newly created quarks (antidiquarks) com¬ 
bines with each of the primary c quarks to respectively 
hadronize into two open charmed hadrons, which can be 
described by an assigned concrete hadronization model 
(for details, see 19413 ). We adopted the Lund String 
Model (LSM) [^, Il0l| in our previous work . In this pa- 
per , we employ the Quark Combination Model (QCM) 
|13l Il4| to deal with the hadronization and compare with 
the corresponding previous results in Q to study to what 
extent different hadronization models interplay with this 
kind of color connection. As other hadronization mod¬ 
els, QCM successes in reproducing experimental data of 
e+e“ —5> h's and pp{p) -T h's processes. As a matter of 
fact, the baryon to meson ratio [I10, and constituent 
quark number scaling of elliptic flow V 2 measured at 
RHIC experiments can be naturally understood by QCM. 
Hence, such comparisons of various hadronization models 
connected with the same color connection can shed light 
on different aspects of soft strong interactions reflected 
by various models. 

The details of the application of QCM are similar to 
that of LSM in Ref. As shown in Fig. [1] e.g., for 5cc 
production, the distribution of Sec is described by the 
Peterson formula [l^ 
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where eg is a free parameter which is expected to scale 
between flavors as eg oc 1 /toq. In the following, we 
will show the results corresponding to eg = 1/25 for the 
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FIG. 1: The procedure of Hcc production for the special color 
connection of Eq. ©• Solid circles represent quarks, while 
hollow circles represent antiquarks. Each of two newly created 
quarks and one primary c are connected to form the color 
singlet system in (3). 

(cc) diquark. The complementary antiquark transits to 
an antibaryon by combining with an antidiquark and the 
distribution is described with the help of a standard frag¬ 
mentation function [l^ as well as the analytical formu¬ 
lation employed in PYTHIA [l^ 

f{z) oc — z)°'exp{—hm\/z), (3) 


[2l|). QCM treats meson and baryon formation uni¬ 
formly, so the production ratio of baryon to meson can 
be directly obtained as 

= 3(aiV + 6) ’ 

where a = 0.66 and b = 0.56. 

As in other phenomenological models, some inputs 
must be given to determine the momentum distribution 
of primary hadrons. In this paper, we simply adopt the 
widely used LPSA which is equivalent to the constant 
distribution of rapidity. Hence a primary hadron i is uni¬ 
formly distributed on the rapidity axis, and its rapidity 
can be written as 

Y = Z + 0 < < 1, (6) 

where is a random number; Z and Y are two argu¬ 
ments, and can be determined by energy-momentum con¬ 
servation in such a color singlet system 

H H 

= W, = 0, (7) 

i=l i=l 


where a and b are free parameters. In our program, 
for simplicity we use Eq. @ w ith a = 0.3 GeV“^ and 
h — 0.58 GeV“^ as in Ref. The balancing diquark 
is broken by the interactions within the remaining sys¬ 
tem, then each of two quarks becomes connected to the 
primary c quark to form two q'c separated systems, and 
each is described by QCM [3111 for the hadronization. 

In the framework of QCM, Quark Production Rule 
(QPR) and Quark Combination Rule (QCR) are adopted 
to describe the hadronization in a color singlet system, a 
simple Longitudinal Phase Space Approximation (LPSA) 
is used to obtain the momentum distribution for primary 
hadrons in their own system, and then this hadroniza¬ 
tion scheme is extended to the multi-parton states. For 
QPR, in a color singlet system formed by q'c, the average 
number of the quark pairs < N > produced by vacuum 
excitation is given by 


< >= — a — 1, (4) 

where W is the invariant mass of the system, a = Pfh—j, 
/3 is a free parameter, fh is the average mass of newborn 
quarks, and Mqi and Mg are the masses of endpoint quark 
and antiquark. 

For the combination process, all kinds of hadroniza¬ 
tion models satis fy th e near correlation in rapidity more 
or less. In Ref. [^, we have shown that the nearest 
correlation in rapidity is in agreement with the funda¬ 
mental requirements of QCD, and determines QCR com¬ 
pletely. The rule guarantees that the combination of 
quarks across more than two rapidity gaps never emerges 
and that N quarks and N antiquarks are exactly ex¬ 
hausted, and thus the unitarity is guaranteed (see Ref. 


where Ei and P^i denote the energy and the longitudi¬ 
nal momentum of the Ith primary hadron respectively, 
obtained by 


771 exp{Yi)-\-exv{ — Yi) 

E^ = rriTi ^ ^ 

P , exp(Yi)—exp( — Yi) 

Lt = rriTi ^ ^ 

where rriTi is given by 


( 8 ) 


rriTz = 


Pn 


(9) 


where rrii is the mass of the Ith primary hadron, and Pxi 
obeys the distribution 


^ ^ H H ^ 

/(Pti, • ■ •, Pth) oc Wexp{ - 

i=l ^ i=l 

In this paper, we set a = 0.2 GeV. Eq. (HOD is just what 
LSM uses. 

For the numerical results, similarly to [^, the center- 
of-mass energy is set to Z° pole. At first, we investigate 
the hadronization effects on jets. Here we take JADE 
algorithm as j/y = to define jets. The pa¬ 

rameter yc.ut is thus introduced and two partons/particles 
are considered as being in one jet when < ycut- We 
apply the jet algorithm to the parton level to obtain the 
three-jet events of e+e“ —^ (cc)cc. Hereafter, three-jet 
events are defined with y^ut = 10“^ if no explicit expla¬ 
nation is given. Fig. [2] shows the energy fraction distri¬ 
bution of the jet involving and that of the jet from 
the primary c. It is obvious that the jet involving Sgg 
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is the hardest one. Both QCM and LSM give the simi¬ 
lar results. In the hadronization process, the energy and 
momentum are broadened relative to those of the par- 
ton level. To show this effect, the distribution for the 
invariant mass of two primary c jets system is displayed 
in Fig. [3][a). The distribution spreads more widely than 
that of two primary c quarks. 




FIG. 2: The distribution of the energy fraction for the jet 
involving compared to the jet from the primary c as a 
function of the scaled dimensionless variable xe = 2 EI\/s. 
The solid line represents the jet involving Hcc, and the dotted 
line represents the jet from the primary c. (a) / (b) stands for 
QCM/LSM. 




FIG. 3: (a) The distribution for the invariant mass of 

two primary c quarks compared to two primary c jets sys¬ 
tem as a function of the scaled dimensionless variable xm = 
2Mir,vI^/s. The solid black line represents two c quarks, and 
the dashed red (dotted green) line represents two primary c 
jets system for QCM (LSM); (b) The transverse-momentum 
distribution for Hcc compared to the diquark (cc) as a func¬ 
tion of the scaled dimensionless variable xt = 2ptI\/s- The 
solid black line represents (cc), and the dashed blue (dotted 
green) line represents Sec for QCM (LSM). 

To further demonstrate the hadronization effects, the 
transverse-momentum with respect to the thrust axis dis¬ 
tributions of 5cc are shown in Fig. [3Kb). The thrust axis 
is determined from the final hadrons system. The abso¬ 
lute values of the distributions more or less depend on 
the parameters in Eq. m, ® and the entire hadroniza¬ 
tion. Once data are available at the future relevant ex¬ 
periments, e.g., International Linear Collider, the Higgs 


factory and factory, etc., the fragmentation functions 
and the parameters can be tuned according to the com¬ 
parisons with data. From Fig.|3Kb), one can also find that 
there is little difference between the results from QCM 
and LSM. 

In the special color connection discussed in this paper, 
it can be clearly seen from Fig. [T] that there is no color 
flow between two primary c quarks. So after hadroniza¬ 
tion, there should be few hadrons emerging in the phase- 
space between them. This is the so-called ‘string effect’ 
and the corresponding event shape should not be sensi¬ 
tive to which hadronization model is employed. This is 
clearly demonstrated by what is shown in Fig. [d] Here 
we choose the more symmetric three-jet events by requir¬ 
ing that the angle between any two jets is larger than 
7 r/2. Because of momentum conservation, all three jet 
momenta must be in the same plane {V) in the e“''e“ 
center-of-mass frame . The three-momentum of each 
final-state particle ki is projected onto one of the three 
regions between the jets to obtain the two-dimensional 
vector in the plane V. The three-momentum of the jet 
that contains Sec is chosen to be the x axis. The angle 
between and the x axis is the azimuthal angle (j) of 
the corresponding particle. We can then calculate the 
final particle-number (energy) distribution l/NdN/d(j) 
{1 /E dE/d(j)). The corresponding results are shown in 
Fig. m Obviously, both QCM and LSM can give this 
kind of effect. 




FIG. 4: (a) The solid (dashed) line represents for QGM 

(LSM); (b) The solid (dashed) line represents for QCM 

(LSM). The symmetric three-jet events are defined here to 
require each relative angle to be larger than 

The above distributions of the transverse-momentum 
and final particle-number (energy) as well as the string 
effect which is the fingerprint of the fact that there is 
no color flow between the separate primary heavy anti¬ 
quarks, are found determined by the color configuration 
and not sensitive to the hadronization models. Therefore 
the study of color connections is important for under¬ 
standing the non-perturbative quantum chromodynam¬ 
ics. The future relevant experiments, e.g., International 
Linear Collider, the Higgs factory and factory, etc., 
can provide opportunities for this study. Moreover, it is 
well known that the hadronization models should be uni- 
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versal for different hadronization processes. So QCM and 
LSM, etc., can also be applied in other hadronization pro¬ 
cesses, e.g., pp collision in Large Hadron Collider (LHC) 
and DIS. For the case that four-quark system emerges 
which has large rapidity gap with other clusters, the 
hadronization effects can also be studied by QCM when 
suitable observables introduced, e.g., the ratio of baryon 
to meson of Eq. In addition, (cc) can be regarded 
as a colored cluster, so if similar heavy colored parti¬ 
cles [^ . [^ beyond Standard Model, e.g., hidden valley 
quarks |26l. are produced at LHC, our method can 
provide useful hints to investigate the related hadroniza¬ 
tion phenomena. 
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